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Abstract. We have investigated the influence of several 
line overlaps between 28 SiO , 29 SiO and 30 SiO on maser 
emission from the envelopes of late-type stars using the 
Large Velocity Gradient radiative transfer formalism mod- 
ified to account for overlaps. Our model includes: (i) 40 
rotational levels for each of the v = to 4 vibrational 
states and each isotopic species; (ii) collisions with molec- 
ular and atomic hydrogen; and (iii) radiative pumping by 
a shell of dust grains surrounding the SiO masing layers. 

We have shown that several line overlaps may play a 
role on the excitation of the 28 SiO v = 3, J=1^0 maser, 
and on the anomalously weak u = 3, J = 2 — > 1 emission. 
We suggest that 28 SiO maser emission in the v = 4, J = 
5 — > 4 line, the only one observed in this vibrational level, 
may also result from line overlaps. For 29 SiO , infrared 
line overlaps successfully explain observed maser emission 
in the v = 0, J = 5 — > 4 and u = 1, 2, J = 6 — > 5 lines. 
30 SiO , v = 0, J = 2 -> 1, and v = 0, J = 4 -> 3 maser 
emissions could also result from line overlaps, and we find 
that several overlaps could explain the newly discovered 
maser line 30 SiO u = 1, J = 1 — s-0. Finally, we make some 
predictions for new 28 SiO , 29 SiO and 30 SiO maser lines. 
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1. Introduction 

SiO masers are excited in the inner layers of circumstel- 
lar envelopes of late-type stars, before the dust formation 
zone, but beyond the photosphere where outflow and in- 
flow of matter or shocks coexist and give rise to complex 
physical conditions. 

Since the discovery of SiO maser emission by Snyder & 
Buhl (1974), many papers have been published to explain 
the SiO maser phenomenon, and several improvements in- 
volving collisional and radiative pumping have been pro- 
posed to model the emission from late-type stars. How- 
ever, it is still difficult to explain the disparities between 
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some rotational line intensities observed inside a same vi- 
brational state or, more generally, in adjacent transitions 
of 30 SiO , 29 SiO and v > 3 28 SiO . Olofsson et al.(1981, 
1985) introduced the hypothesis of a line overlap between 
two ro-vibrational lines of SiO and water to explain the 
weakness of the 28 SiO v — 2, J = 2 ^ 1 emission. This 
emission was observed by Bujarrabal et al.(1996) in stars 
with different O-abundances, and the overlap between two 
infrared lines of H2O and SiO was confirmed. In addi- 
tion, several overlaps between the infrared lines of the 
main ( 28 SiO ) and rare ( 29 SiO and 30 SiO ) isotopic species 
of silicon monoxyde have been suggested to play a role 
in the excitation of these species. Using the LVG ap- 
proximation Cernicharo et al.(1991) modeled the pump- 
ing of 29 SiO by 28 SiO . The importance of overlaps was 
further investigated by Cernicharo & Bujarrabal (1992) 
and Cernicharo et al.(1993). Recently, Gonzalez- Alfonso et 
al.(1996) and Gonzalez- Alfonso & Cernicharo (1997) used 
a non-local radiative transfer code to study the question 
of line overlaps. They successfully explained most of the 
29 SiO , 30 SiO and high-v (v > 3) 28 SiO masers; however, 
several lines remain unexplained. 

Our main goal in this work is to develop a simple and 
modular calculation tool in order to investigate the im- 
pact of various line overlap effects among the three main 
isotopic forms of SiO. We do not intend here to realis- 
tically model the SiO maser line profiles. We rather are 
interested in predicting and discussing relative peak line 
intensity ratios. In Sect. 2 we give details on spectroscopy, 
the used radiative transfer code and physical parameters, 
and on validation of our code. In Sect. 3, we include several 
line overlaps and compare our results to the observations. 
Some major results of this work are summarized in Sect. 
4. 

2. The Model 

The LVG code used here must be regarded as a simple 
but powerful computational tool for maser emission pre- 
diction. Bujarrabal (1994) has clearly demonstrated that 
the LVG code applied to the SiO maser zone is a fairly 
good approximation in comparison with the exact radia- 
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tive transfer solution. Other uncertainties do not play a 
minor role in comparison with the differences between the 
LVG and the exact radiative transfer results; in particular 
uncertainties in the collision laws and the exact geometry 
of the maser region result in weaker predictions. Thus, we 
believe that the physical laws and the spherical symmetry 
adopted in our model are sufficient for the LVG approxi- 
mation. 

2.1. Spectroscopy, collisions and radiative transfer 

Modeling the SiO emission requires the use of accurate 
spectroscopic data and collisional rates. Energy levels of 
SiO have been calculated using the most recent calcula- 
tions of the Dunham coefficients (Campbell et al. 1995). 
For the Einstein A-coefficients we used the recent work of 
Drira et al.(1997). 

Basic data concerning the rotational and rovibrational 
rate coefficients of SiO with H 2 and their dependence on 
temperature were taken from Bieniek & Green (1983). We 
also used the new rate coefficients derived by Lockett & 
Elitzur (1992) for v > 3 and Av = 2 transitions and we ap- 
plied corrections to their rates to compensate for the effect 
of the limited number of rotational and vibrational levels 
taken in the calculations. In addition, Langer & Watson 
(1984) made a first attempt to estimate collision rates of 
SiO with atomic hydrogen from a comparison with the 
expected H-CO and H 2 -CO rates. The main features of 
these H-SiO rates used in our work arc: the vibrational 
collision rates are 10 times larger than those for H2, and 
the rotational collision rates are decreased by a factor of 
7. H can be formed in two different ways: (i) the stellar 
radiation may photodissociate H2 close to the star, thus 
forming an atomic hydrogen layer; (ii) shocks may dissoci- 
ate Hi to produce a uniform distribution of H throughout 
the envelope. We assume here that x[Hi/H2\ is constant 
and equal to 10~ 2 as in Bowers & Knapp (1987). Adding 
atomic hydrogen to H 2 increases the role of collisions. Note 
however that its impact on the inversions is weak (a few % 
for x[H/H 2 ] < 1CP 1 ). Beyond this value, inversions for the 
lowest rotational levels tend to decrease contrary to what 
is observed for high J transitions. The higher the vibra- 
tional level, the larger the effect of the atomic hydrogen 
is; this is expected because of larger vibrational collision 
rates for H than for H2. Nevertheless, we stress that these 
collision rates and the relative abundance of the atomic 
hydrogen are poorly known. 

Because the purpose of this work is to explain in a 
simple way the main characteristics of 28 SiO , 29 SiO and 
30 SiO emissions, we have used the LVG (Sobolev 1958) 
approximation, based on the model of a homologously ex- 
panding envelope, in which the escape probability does 
not depend on the angle. If the velocity gradient (loga- 
rithmic velocity gradient e r = d In V/d In r) in the circum- 
stellar envelope is large enough, each cell of the discretized 
medium will be decoupled from all other cells, thus allow- 



ing a local treatment of the radiative transfer. The in- 
teraction area is limited to a small zone where photons 
can be absorbed or emitted around a resonance point. We 
have limited the maximum size of the cells to the Sobolev 
length, L, defined by 

AV D /V 
L = r 

£r 

with 

AV D = \j AVth 2 + AV turb 2 

where r is the radial distance, AV t h the thermal veloc- 
ity, AVturb the turbulence, and V the expansion velocity. 
AVd characterizes the local absorption coefficient. The 
formula giving the Sobolev length applies to any form of 
expansion velocity and to the general case in which the 
photon escape probability depends on the direction. For 
a spherically symmetric envelope and for the special case 
where V is proportional to r, the angle-dependent terms 
vanish, the logarithmic gradient e r = 1, and the escape 
probability is isotropic and proportional to the inverse of 
the peak opacity (as soon as the opacity becomes large). 
The Sobolev length is just the length over which the ve- 
locity varies by an amount corresponding to the width of 
the local absorption coefficient. One essential advantage 
of the LVG code is that it requires less computer time 
than more exact treatments of the radiative transfer prob- 
lem based on integral solutions (Bujarrabal 1994) or on 
the Monte Carlo method (Gonzalez- Alfonso & Cernicharo 
1997). In the frame of the LVG formalism and for phys- 
ical conditions relevant to the circumstcllar environment 
of late-type stars we have derived the solutions of the sta- 
tistical equilibrium equations and the opacities of various 
SiO lines. More complex effects related to non spherical 
geometry, polarization or saturation and beaming angles 
are beyond the scope of this paper. Concerning satura- 
tion, we emphasize that our rate equations are solved for 
both positive and negative opacities, but that beaming re- 
lated to saturation cannot be investigated here because we 
assume the escape probability to be isotropic. Note that 
beaming angles for saturated masers in the case of the 
LVG approximation remain an unsolved problem (Elitzur 
1992). Saturation and the effect of competitive gains in 
SiO lines were investigated in Doel et al.(1995). 

Our calculations include 40 rotational levels for each 
of the 5 vibrational states v — to 4. To correct for the 
limited number of rotational and vibrational levels (this 
overestimates the opacities), we have applied a correction 
similar to that used by Bujarrabal & Ricu (1981) for the 
rotational populations. The required final precision on the 
populations is better than 10~ 4 . Calculations are done 
simultaneously for 28 SiO , 29 SiO and 30 SiO . We derive the 
populations for each level (v, J) and derive the opacities 
for all allowed transitions. 
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2. 2. Circumstellar parameters 

We have used physical parameters appropriate to the cir- 
cumstellar environment of evolved stars (see details in 
Herpin 1998). 



* Kinetic temperature 

The gas temperature follows the law adopted by Langcr 
& Watson (1984): 



T 



0.5 



where i?* is the stellar radius (R± = 7.7 10 cm), T* is 
the stellar temperature (T* = 2500 K) and r is the radial 
distance. This law is in rough agreement with the obser- 
vations of SiO thermal emission of Lucas et al.(1992) and 
Bujarrabal et al.(1986), and also with the calculations of 
Willson (1987) who finds 1500 K in the maser region, for 
a typical Mira. Turbulence is a free parameter in our code 
and contributes to net line broadening. 



This abundance decreases with the distance to the star, 
the SiO molecules being progressively condensed onto 
the grains. We follow the law adopted by Bujarrabal et 
al.(1989), and we adopt for the 29 SiO and 30 SiO isotopic 
abundances (Ai) relative to 28 SiO , 1/19.5 and 1/29.5, re- 
spectively (see fig.0): 

X(r) = Xo[l - ^ + X'UA 

with xo the initial ratio (here the solar abundance), and 
x' the remaining fraction of SiO after the grain formation 
(10 -2 ). The terms n(r) and kq depend respectively on the 
radiative pressure and the total quantity of grains formed. 
We take: 



n(r) = kq 



{r-Rif 



{r-Rxf + {Rd-Rif 



where Rd is the radius where half the total quantity of 
grains is formed (1.5 10 15 cm), and Ri is the initial point 
(2. 10 14 cm). 



* Expansion velocity 

The expansion velocity of the gas adopted in our model 
is: 



V = V X 



1 



1 + ( i^ )e „ 



where R acc = 1.8 10 15 cm, e v £ [1; 3] and is the asymp- 
totic expansion velocity of the gas (9.5 kms -1 ). A logar- 
rithmic velocity gradient of 3 was used in the calculations. 



* ^density 

We adopt (Elitzur 1992): 



M 



n,H 2 = 



8-7T r 2 Vm v 



with V the expansion velocity, M the steallar mass loss 
rate and m p the proton mass (see FigJ|). Maser effects 
will require densities of a few 10 9 particles/cm 3 , whereas 
thermalization is clearly reached beyond 10 10 — 10 11 
particles/cm 3 . 



X(SiO/H 2 . 




r(cm) 



Fig. 1. SiO/H2 abundance x (continuous line) and H2 
density (in cm -3 , dotted line) versus the distance to the 
star r(cm) for X o = 5 10~ 5 , £jf = 5.3 10~ 7 M /yr and 
R* = 7.7 10 13 cm. 



* SiO abundance 

The SiO/Bi2 abundance ratio is uncertain, and we adopt 
the solar abundance by default (xo[SiO/H 2 ] = 5. 10~ 5 ). 



cmj 

A \-p 



* Grain radiation 

IR spectra from the envelopes of late-type stars show that 
circumstellar grains have various chemical compositions. 
Here we only consider silicate grains and we assume that 
most circumstellar grains lie in a shell outside the SiO 
maser zone (with Ri n = 10 15 cm, and R ou t = 10 16 
The extinction coefficient of the dust is Q v ~ Qo{j-) 
with p G [1;2] for wavelengths A > lfxm (Mathis 1990). 
For silicate grains, we take p = 1.1, Ao = 80/xm, and 
Qo = 2.10 -3 according to the work of Ivezic & Elitzur 
(1995). We consider that the silicate grains behave as 
black-bodies with temperature Td- We assume that the 
dust is optically thin, and that Td is constant throughout 
the envelope (Td =600 K). We take, following Netzer & 
Elitzur (1993), Xd = = 1CT 2 , 3.0 g.cm~ 3 for the 

grain volumic density, and 5.10 _6 cm for the grain radius. 

2.3. Line Overlaps 

In this work, we only deal with local line overlaps due to 
thermal line broadening and local turbulence in successive 
circumstellar layers. For simplicity, two nearby spectral 
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lines are considered to overlap when AV < AVo where A V 
is the difference in velocity between the line centers, and 
AVo is an a priori fixed value. In the examples discussed 
in Sect. 3 we have adopted AVo = 5 km s~ x corresponding 
to turbulent inner layers (AVthermai — 1-2 kms^ 1 ). Line 
overlaps are treated in a very simple manner. For two 
nearby lines v\ 2 and ^34, we derive the opacities T12, T34, 
and the individual source functions S12, S34,. Adding the 
absorption and emissivity coefficients of each line, the 
source function becomes 



and the total opacity is now t\2 — r i2 + T34- The average 
intensity is then given by: 

Ju = [1 - /?12 V)]^' + 0i2bb'(r')I bb 

+{3l2bg'(T')(I bg + Idust) 

where is the usual escape probability. The cosmic back- 
ground h g is described by a 3 K blackbody and the dust 
contribution is given by Idust ■ The emission of the central 
star Ibb is described by a blackbody spectrum at temper- 
ature T*; a geometrical dilution factor is applied. In these 
formulae, the populations of the 4 levels are involved. If 
T12 < T34, after overlap we have S\2 = S34 ~ S34 and 

~ T34; this implies minor changes for T34 and the 77.3 
and m populations. On the other hand, photons from the 
stronger line will be absorbed by the more optically thin 
line 1 — 2, thus enhancing the non-equilibrium distribution 
of the populations in the levels 1 and 2. 

Such a treatment of the overlaps is oversimplified, but 
allows us to test quickly the effects of this mechanism. 
Our code incorporates three overlaps simultaneously; this 
number can of course be increased. To account for line 
overlaps, the code first derives the populations and opaci- 
ties separately for each isotope without overlap, and then 
solves again for the radiative transfer, combining the over- 
lapping lines. We describe in Sect. 3 how we proceed to 
test the effects of line overlaps. 

In fact, non local line overlaps have long been recog- 
nized as important in the excitation of the OH radical, 
and, more recently, in SiO models (cf. Introduction and 
Sect. ^). Line overlaps may occur even if two lines are not 
excited in immediately adjacent gas layers; e.g. a Doppler- 
shifted line emitted from a first cloud may interact with 
another line in a second cloud provided that the two clouds 
have the appropriate relative velocities. In a somewhat dif- 
ferent process, Olofsson et al.(1981, 1985) suggested that 
the anomalously weak v = 2, J = 2 — > 1 emission line 
in stars could be explained by non local overlap with the 
H 2 V2 1275 — > v\ 11 6 6 line. Photons emitted by the wa- 
ter line increase the radiation field at the frequency of the 
vibrational transition of SiO creating an excess of absorp- 
tion for this SiO line, and thus destroying the inversion in 
v = 2, J = 2 — ► 1. This mechanism was recently investi- 
gated in detail by Bujarrabal et al.(1996). 
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Fig. 2. 28 SiO negative opacities for the transition J = 
1 -> in v = to 4. We have taken T* = 2500 K, 
R+ = 7.7 10 13 cm, the logarithmic velocity gradient e = 3 
and ^ = 5.3 10~ 8 M Q /yr. (The H 2 density is fixed to 10 9 
cm -3 at a distance of 2.8-/?*.) V is the expansion velocity. 

2.4. Model validation 

To test the validity of our code, we have used the same 
initial conditions as in Alcolea et al.(1989). (For this com- 
parison, we have not applied corrections on the collision 
rates for the limited number of levels.) The main differ- 
ences with Alcolea et al. are the number of rotational 
levels (40 here compared to 12 in Alcolea et al.), the colli- 
sion rates (Langer & Watson 1984 in Alcolea et al.), and 
the more recent spectroscopic coefficients used here. In 
Figure ^, we plot the J = 1 — > negative opacities for 
v = to 4 versus the ratio n(SiO)LV~ 1 e~ 1 (where L is 
the Sobolev length, V the expansion velocity, e the loga- 
rithmic velocity gradient). Our results are much similar to 
those of Alcolea et al. and show that the J = 1 — > inver- 
sions are well differentiated according to the vibrational 
state. However, some minor discrepancies remain. In par- 
ticular, our work shows that v = 1 inversions occur over a 
broader range of parameters. We do not obtain increasing 
opacities with increasing v; but in Alcolea et al. this trend 
was probably the consequence of their smaller number of 
rotational levels. Finally, we note that results in Fig. 2 are 
similar to those of Lockett & Elitzur (1992). 

3. Discussion of line overlap effects 

Standard radiative and collisional SiO pumping models 
are able to predict several of the observed features of 
SiO maser emission from late-type stars (e.g., Bujarra- 
bal 1994). These models predict similar line profiles and 
smooth line strength variations when one goes up the ro- 
tational energy ladder. They cannot explain, however, the 
anomalous line intensity ratios observed among several ad- 
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Fig. 3. Negative opacities of newly explained maser lines for 29 SiO (top) and 30 SiO (bottom) with (continuous line) 
and without (dashed line) overlap, versus the distance to the star (one stellar radius = 7.7 10 13 cm). The overlaps used 
in these computations are given in bold face in Table [l]. The H% density and the abundance x[^iO / H%\ are respectively 
10 9 cm~ 3 and 5 10~ 5 at a distance of 2.8i?*. 

jacent rotational transitions in high v states of 28 SiO and al.1996, Gonzalez- Alfonso & Cernicharo 1997) includes a 
in v — to 3 states of 29 SiO and 30 SiO . On the other non-local treatment of line overlap effects. Their model 
hand, line overlap effects may strongly modify the SiO correctly predicts several observed features. It also pre- 
populations and thus explain apparent line anomalies. The diets excitation of the 29 SiO v = 3, J = 8 — > 7 transi- 
most complete work in this domain (Gonzalez- Alfonso et tion which was subsequently detected (Gonzalez- Alfonso 
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Table 1. Possible explanations of observed maser lines not explained in previous works. The overlaps used in Figs. 
U and [I| are given in bold face. Note the 30 SiO v = X, J =1 — > maser line newly discovered by Cho & Ukita (1998). 



Explained maser lines 


Overlaps 


2y SiO v = 0, J = 


5 


— > 4 (weak) 


. 2y SiO v = 


1-0,J = 


1 - + 2S SiO « = 2 


— 1, J = 4 — 


3 


2y SiO v = 1, J = 


6 


-> 5 


. 28 SiO ti = 


2-l,J = 


4 - 3 + 2H SiO « = 1 


- 0, J = 1 - 











and 28 SiO v 


= 2 — 1, J 


= 5 — 4 + SiO i! = 


1 — 0, J = 2 


— 1 








. 28 SiO u = 


3 - 2, J = 


2 - 3 + 29 SiO u = 2 


-1,J = 5- 


6 








. 29 SiO v = 


= 2-1,3 


= 15 - 14 + 29 SiO 


v = 1 - 0, J 


= 6-5 


2y SiO v = 2, J = 


6 


-> 5 


. 3U SiO v = 


= 1 0, J 


= 1 - + 2y SiO v = 


= 1 0, J = 


3 4 








. 28 SiO v = 


3-2, J = 


2 - 3 + 29 SiO v = 2 


-1,J = 5- 


6 


au SiO v = 0, J = 


2 


-» 1 


• 3U SiO v = 


= 1 0, J 


= 1 - + 29 SiO v = 


= 1 0, J = 


3 4 








. 28 SiO v = 


2-l,J = 


3 - 4 + 30 SiO w = 1 


— 0, J = 1 - 


2 


' iu SiO « = 0, J = 


4 


-> 3 


• ^SiO v = 


= 2 1, J 


= 3 - 4 + ;iu SiO v = 


= 1 0, J = 


1-2 


■ lu SiO?; = 1, J = 


1 


-> 


. 3U SiO v = 


1-0,J = 


1 - + 2a SiO v = 3 


-2, J = 11 - 


- 10 








and 29 SiO v 


= 3-2, J 


= 11 - 10 + 28 SiO u 


= 2 - 1, J = 


= 2-3 








. 29 SiO v = 


= 1 0, J 


= 3 4 + 30 SiO v = 


= 1 0, J = 


1-0 








. 28 SiO v = 


2 - 1, J = 


3 - 4 + 30 SiO v = 1 


- 0, J = 1 - 


2 



& Cernicharo 1997) thus demonstrating the importance 
of line overlaps. Nevertheless, excitation of several transi- 
tions remains unexplained in the work of Gonzalez- Alfonso 
et al. In particular, emission from 29 SiO v = 0, J = 5 ^ 4 
and v = 1,2, J = 6 -» 5 and 30 SiO v = 0, J = 2 -> 1 and 
4^3 must be understood (see discussion below and Table 
0)- 

The present model is conceived as a simple tool to 
quickly investigate the impact of various line overlaps be- 
tween ro- vibrational lines of 28 SiO , 29 SiO and 30 SiO . The 
model incorporates up to 3 simultaneous different over- 
laps, namely six lines, ac cording to the simple prescrip- 
tion outlined in Sect. 2.3. Once the maximum line shift 



allowed betwee n tw o lines is fixed (this depends on AVb 
defined in Sect. 2.3), a specific module of our code searches 
for all possible overlaps. Within the 10 first rotational 
transitions we find a total of 27, 48 and 79 overlaps for 
AVb = 5,10 and 15 kms^ 1 , respectively. There are of 
course more overlaps involving higher J rotational lev- 
els, and we have also explored the influence of high J- 
level transitions overlapping low J-level transitions. As 
briefly explained in Sect. 2.3, the effect of a line over- 



lap between two lines strongly depends on the relative 
strength of these two lines. In general, when an overlap 
occurs between 28 SiO and 29 SiO or 30 SiO transitions, the 
rare isotopic transitions which are optically thin tend to 
be more affected than the 28 SiO transitions; and the im- 
pact on the more optically thin lines is higher for lower 
v states. Our code derives the net optical depths com- 
pared to the non-overlapping cases. Taking AVb = 5 
km s _1 results in 31 overlapping line pairs. These overlaps 
affect 28 SiO , 29 SiO or 30 SiO rotational transitions with 
v < 4 and J < 10. There is however the interesting case 
of 28 SiO u = 2-l,J = 3->4 which requires AVb - 10 
kms _1 to overlap with 30 SiO u = l-0, J=l->2. 

Because of radiative pumping through high v levels, we 
stress that very different (v, J) levels are connected. That 



is why some infrared overlaps have sometimes a strong in- 
fluence on ro- vibrational levels not involved in these over- 
laps: because of this excitation mechanism numerous very 
different levels are connected together. 

We have grouped our results into four categories: (i) 
the overlaps invoked by Gonzalez- Alfonso et al.(1996) and 
Gonzalez-Alfonso & Cernicharo (1997) to explain some 
important newly discovered rotational transitions; (ii) the 
overlaps which may account for the behaviour of maser ro- 
tational transitions explained elsewhere by other models 
and other overlaps; (Hi) the overlaps providing explana- 
tions for the maser lines not yet explained; and (iv) the 
overlaps suggesting a search for new maser lines. 

• In the first category it is worth mentioning that the 
overlap of 28 SiO v = 2 - 1, J = 20 - 21 with 29 SiO v = 
3 — 2, J = 8 — 9, although it involves relatively high v 
states for 29 SiO , results in strong enhancement of the 
29 SiO u = 3,J = 8 — 7 transition. This was also predicted 
by Gonzalez- Alfonso et al.(1996). We have also verified 
that two line overlaps, 28 SiO v = 2 — 1, J = 2 — 3 with 
29 SiO v = 3 - 2, J = 11 - 10, and 29 SiO v = 3 - 2, J = 
11 - 10 with 30 SiO v = 1 - 0, J = - 1, result in a 
clear enhancement of the 29 SiO v = 3, J = 11 — ► 10 
line intensity. Another example is the enhancement of 
the 30 SiO u = 2,J = 8^7 transition resulting, as in 
Gonzalez- Alfonso et al.(1996), from the effects of 3 simul- 
taneous line overlaps. 

• We have successfully explained (case (ii)) the en- 
hancement of the 28 SiO w = 3,J=l— >0 and v = 4, J = 
5^4 transitions, and of the 30 SiO v = 0, J = 5^4 
transition. Standard modelisation does not easily invert 
the 28 SiO u = 3, J=1^0 line although it is ob- 
served in a variety of sources (Scalise & Lepine 1978, 
Alcolea et al.1989). One overlap may lead to the exci- 
tation of this line: 28 SiO v = 3 - 2, J = 2-3 with 
29 SiO v = 2 — 1, J = 5 — 6. On the other hand, obser- 
vations indicate that the 28 SiO v = 3, J = 2 — 1 line 
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Table 2. Predicted new maser transitions. The lines in 
bold face have not been searched for maser emission as far 
as we are aware. 



2H SiO 


V 


= 4 


J = 


:i -» 


2 








J = 


4 -> 


3 


M SiO 


V 


= 


J = 


4 -> 


3 




V 


= 1 


J = 


2 -» 


1 








J = 


5 -> 


4 




V 


= 2 


J 


1 -> 











J = 


3 -> 


2 




V 


= 3 


J = 


2 -» 


1 








J = 


4 -> 


3 


:lu SiO 


V 


= 1 


J = 


2 -> 


1 




V 


= 2 


J 


2 -> 
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should be quenched. This could result from four differ- 
ent overlaps including the low J rotation level overlap of 
28 SiOv = 2-l, J = 3-4 with 30 SiOu = 1-0, J = 1-2. 
This overlap is also consistent with the extinction of the 
u = 3, J = 2 — 1 line and with the enhancement of the 
v — 4, J = 5 — > 4 line as observed by Cernicharo et 
al.(1993) in VY CMa. Our model easily leads to maser 
emission of the v — 0, J = 1 ^ transition. This is ob- 
served in several stars for 29 SiO (Alcolea & Bujarrabal 
1992) and 30 SiO (Cernicharo & Bujarrabal 1992). 

• Table |j] (case (ra)) and Fig.|| suggest explanations of 
the 29 SiO and 30 SiO maser transitions not explained in the 
work of Gonzalez- Alfonso et al.(1996). Only overlaps can 
explain inversions of the 30 SiO v = 0, J = 2 ^ 1,4 — ► 3 
(not explained in previous works) or J = 5 — > 4 tran- 
sitions. Line overlaps are also required to explain sev- 
eral 29 SiO v = 1 and 2 rotational transitions. For the 
v = X, J = 6 — ► 5 line observed by Cernicharo & Bu- 
jarrabal (1992), Table [l] suggests that it may be ex- 
plained by several overlaps. Of special interest is the si- 
multaneous overlap of 28 SiO v = 2 — 1, J = 4 — 3 with 
29 SiO v = 1-0, J= 1-0 and of 28 SiO v = 2-1, J =5-4 
with 29 SiO v — 1 — 0, J = 2 — 1. The same arguments ap- 
ply to v — 2, J = 6 ^ 5 which has also been observed 
by Cernicharo & Bujarrabal. The overlap of 29 SiO v = 
1 - 0, J = 3 - 4 with 30 SiO u = l-0,J=l-0is 
clearly important because it enhances simultaneously the 
observed J = 4 — 3 and 2 — 1 transitions. It is interest- 
ing to note that we have succeeded in producing weak 
maser emission in 29 SiO v — 0, J = 5 ^ 4 (never ex- 
plained before). Another important result concerns the 
30 SiO v = X, J = X — ► maser emission produced by our 
model with different overlaps (see Table [l] and Fig. ^); this 
line was recently discovered by Cho & Ukita (1998) in TX 
Cam. 

• Table |^ lists some relatively strong maser transi- 
tions as predicted by our model (case (iv)). Note that 
the 29 SiO and 30 SiO u = l,J = 2— >1 maser lines have 
also been predicted by Rausch et al.(1996) on the basis 
of an optical pumping model. Four transitions ( 28 SiO v = 



30 SiO v=l, J=l-() 




2 2,5 3 3,5 4 4,5 

Stellar radius unit 



Fig. 4. Negative opacities for the 30 SiO u = 1, J = X — »0 
maser line with (continuous line) and without (dashed 
line) overlap, versus the distance to the star (one stellar 
radius = 7.7 10 13 cm). The overlap used in this computa- 
tion is given in bold face in Table [l| The H2 density and 
the abundance \[SiO / H2] are respectively 10 9 cm~ 3 and 
5 10~ 5 at a distance of 2.8i?*. 

4, J = 3 -» 2, 29 SiO v = 1, J = 2 -> X and v = 3, J = 
4^3, 30 SiO v = 2, J = 2 — > 1) are easily excited even 
without line overlaps; however, all other transitions re- 
quire line overlaps for being detectable. For 29 SiO v = X 
we find that 6 different overlaps of ro-vibrational lines 
of 28 SiO with 29 SiO tend to enhance the J = 5 -> 4 
transition. Six transitions (given in bold face in Tabic |[ 
28 SiO v = 4, J = 3 -> 2, 29 SiO w = 2,J=1^0,3^2 
and v = 3, J = 2 -> 1, 4 -> 3, and 30 SiO v = 2, J = 2 -> 1) 
have not yet been searched for maser emission as far as 
we are aware. However, all other transitions in Table 2 
have already been observed without success by several au- 
thors (e.g., Cernicharo & Bujarrabal 1992). Nevertheless, 
further investigations of these lines would be useful, espe- 
cially as 29 SiO u = l,J = 2— »1 was recently detected in 
one star (J.R. Pardo, private communication). 

Finally, we stress that other types of line overlaps may 
simultaneously play a role in the anomalous excitation 
of 28 SiO , 29 SiO and 30 SiO . The infrared lines of water 
are especially important in this context because water is 
an abundant species in the late-type stars where SiO is 
present. Apart from the B^O-SiO line overlap discussed 
by Olofsson et al.(1981, 1985) and Bujarrabal et al.(1996), 
there are several other infrared transitions of ortho- or 
para-H20 close to ro-vibrational transitions of SiO and 
isotopes which would deserve modelisation. 

4. Conclusion 

We have presented results for a new and simple model 
of SiO masers which incorporates both radiative and col- 
lisional pumpings as well as line overlap effects. Radia- 
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tive pumping includes stellar radiation and circumstel- 
lar grains although their extinction coefficient is not well 
known. Collisional pumping involves molecular hydrogen 
as well as atomic hydrogen (under the form of simply 
guessed H-SiO collision rates). Our model involves 40 rota- 
tional levels for each of the first five vibrational states and 
derives the populations and opacities for the three species 

28 SiO , 29 SiO and 30 SiO , and for a variety of circumstellar 
conditions. Several low and high J rotational transitions 
are easily inverted in the vibrational states v — 1,2,3 and 
4 as observed in many late-type stars. Prediction of rel- 
ative line intensities must account for line overlap effects 
which have been shown to play an important role in the 
excitation of SiO masers. 

Our model incorporates in a very simple way the ef- 
fects of line overlaps between ro-vibrational lines of 28 SiO , 

29 SiO and 30 SiO . Three different overlaps may be treated 
simultaneously, and we have examined the impact on the 
level populations and opacities of 31 overlapping line pairs. 
We have found that line overlaps are important to enhance 
28 SiO v = 3, J = 1 ^ and v = 4, J = 5 ^ 4 maser emis- 
sion or to turn offu = 3,J = 2^1 line emission. Various 
overlaps may lead to weak 29 SiO v — 0, J = 5 ^ 4 line 
emission and strong 29 SiO v = 1 and 2, J = 6 — > 5 emis- 
sion; these lines had not been explained in other works. 
More generally, several 29 SiO v = 1 and 2 maser lines are 
well explained with line overlaps. We have also suggested 
possible explanations to several 30 SiO maser lines includ- 
ing the newly discovered transition v = 1,J = 1 — » 0. 
In addition, we have predicted that several new maser 
transitions for 28 SiO (v = 4), 29 SiO (v = to 3) and 
30 SiO (v = 1 and 2) could result from line overlap effects. 
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